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L. TNTRODUSITILON 


In 1973 Leo and Boncal [1] converted the Oceanics Inc. 6 
DOF Loads and Motions digital computar simulation program 
for the 100-В Surface Effect Ship to reflect the physical 
Enaracteristics of the XR-3 craft. "Since that time there 
has been a continuous effort at the NPS tò update and 


inprove on this simulatioa (2,3,4,51. 


One of the major icawbacks of that program which 
precludes regular use of it for sea state operations 
Simulation is the excessive computational and turn-around 
tines required [6]. As it is presently ldesired t> use the 6 
DOF simulation for design and evaluation studies of the 
proposed 3K Ton SES control System in varying sea state 
conditions, it is impoctant to significantly reduce the 


Simulation's computational time. 


One area where the possibility exists to achieve a major 
reduction in computational time is tus manner in which 
Sidewall forces and moments агг calculated. The present 
form of the 6 DOF progran uses a large number of discrete 
sections (28 per sidawall) with an off-line table look-up 
scheme to determine the forces genearatei іп each discrete 
section. The simple model derived by McIntyre [7] uses two 
sections per sidewall with  verv simplified geonetry that 
permits on-line computation of tha forces and monents. The 
result of aventual incorporation of this sidewall geometry 
and forces calculation into the 5 DOF program should be a 


Significant reduction of computation and turn-around time. 


Another area of great interest in SES Simulation is 
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gaining insight into the dynamic behavior of the craft's 
vertical motion of pitch and heave due to pressire induced 
lift forces generated in the plenum and seals. As stated in 


ref [8] as recently as August, 1978, 


ıtThe overall precision ЗЕ th2 analysis is satisfactory 
for predicting stability and maneuverability. However, 
several areas present very complex or difficult physical 
phenomena to analyze and further effort should be 


invested to improve understanding of these.!! 
The first area for continuing investigation recoaneanded was; 


ae petter determination of pitch-heave dynamic 


characteristics. !! 


An improved knowledge of the dynamic characteristics 
brougnt about by the pressure induced forces and moments 
could lead to a reduction of the computation time 
required for craft simulation. Thasa pressure affects 
introduce a stiffening in the integration of sea state 
induced oscillations and а better understanding of 
processes could lead t» the removal of this stiffness 


and thus reduced CPU time. 


In this Work both of these areas are considered. First 
of all the various craft physical characteristics are 
modeled, resulting in changes ani improvements ta the 
initial efforts by McIntyre [7]. These equations are 
linearized about a steady-state operating point using 
the Taylor Series expansion. The linearized equations 
ara cast into a State Space format and simulatel on the 
araıtal computer. Comparison is made between the 
tine-domain results of the linear model and tha 6 DOF 
model. The State 52322 representation is used to 


generate a signal flow graph. Mason's Gain Rule is 
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applied and the S-domain PODES ко: the craft's 
characteristic equations are located and it is shown 
that some simplification and reduction of system order 
is possible. The resultant S-domain polynomials are 
used to generate Bode Plots to graphically display the 
craft's simplified frequency response. 


Conclusions are drawn and recommendations are made 
concerning the validty of tha model and areas for future 


study. 


Several procedural notes are in order. Because all 
equations are to be solved digitally using the FORTRAN 
computer language and because digital computer text 
processing facilities were used to generate this 
document, all formulas will be presented in FORTRAN 
format. In the thesis gody reference to the ''noniinear 
nodel'*, ''six degres-of -freedou model'' or '*'6 DOF!" 
model all mean the Loads and Motions program as adapted 
to the XR-3 test craft at the Naval Postgraduat2 School 
n3] Reference to the linear model refers to the 


linearized model of the same craft as developed herein. 
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ІІ. NONLINEAR XR-3 MODEL HEAVE AND PITCH EQUATIONS 
A. BACKGROUND 
The Original development work on the simplified 


two-degree-of£ -freedoa model was done by McIntyre (1. 
Several weaknesses in the original model were noted by 


McIntyre in that work. 


The major portion of the work presented here was 
involved in close examination of all aspects of the 
original model, comparison of tae modeling techniques 
used therin with other methods and remodeling of several 


important craft physical characteristics. 


One major area of investigation was tne lack of pitch 
damping in the original model, which resulted in the 
addition of an added-nass effect to the pitcn moment 


equations. 


Additional changes involve a revision of the force 
equations for the rear seal, major changes ia the way 
the plenum chamber is modeled (also added to #12 loads 
and motion 6-dof model DJ ). and a more rigorous 


development of the planing force equation of 0J . 


There are several areas of theory concerning CAB SES 
modeling where empirical data is lacking due to 
difficulties with physical measurements. One such area 


inviting closer inspection was dynamic variatioa of the 
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effective canter of pressure for the plenum gauge 
pressure lift effect. 4 possible approach is suggested 


herein and included in the model. 
В. ASSUMPTIONS 
The following assumptions are made for this constant 


speed, pitch-heave model: 


1. Planing forces are concentrated at the centroid of 


the keel longitudinal cross section. 


2. Total seal force is lump2i into a pressure 
differential acting on the seal face area in contact 


with the water passing beneath the seal. 


3% A2rodynamic forces on tha forward seal are 
disregarded. 
4. Dynamic variations in moment lever arms are 


disregarded for all monents except seals and tha moment 


due to plenum guage pressure. 


5. Dynamic variations in effective plenum roof area 
Will be included. 


6. Any lift effects lue to aerodynamic forces аге 


lumped in with the planing forces (for reasons to ba 
developed in Zhapt IV.) The moments are included in 
plenum gauge pressure moments . 


7. Constant seal leakage area is assumed. 


8. Since this лойгі is for constant surge velocity, 
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drag and thrust forces are considared to be in 


equilibrium and in balance for moments about th2 Y axis. 


С. COORDINATE SYSTEM 


(see fig. 1) 


1. The origin of th2 cartesian coordinate system is 
located at the calm-water line, longitudinally and 


transversly coincident with the craft center of jravity. 


2. The X axis is the longitudinal coordinate defined 


as positive as X increases towards the bow of tae craft. 


2. The Y axis is the transverse coordinate increasing 


positively towards the starboard side. 


ц. The Z axis is the vertical axis, iacreasing 


positively in a downwari direction. 


5. ZS is the vertical distance from the ОЕ 


{ } 
(D 
2 
сї 
(D 
tj 


gravity to the keel. 


” 


6. Z2 is the vertical distance from the cala water line 


to the vertical center 2£ gravity. 


7. THETA is the pitch angle defined as positive as the 
bow pitches up. 


8. 2ositive pitch moments ar= 1е:1л24 аз попза 5 which 


сапа to cause a positive increase in the angle THETA. 


D. SIDEWALL MODELING 
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As it was desired to keep all physical dimensions 
consistant with the six degree of freedom mol21, scale 
model drawings were dev2loped with the data from the 
sidewall subroutine ОЕ that model. All sidewall 
dimensions usd in this nodel wera taken from those 


drawings. 


The sidewalls are modeled in two Sections, one forward 
of the Wong tudinal 5.6. and the other aft. Fach 
sidewall is considered t> be of constant width at the 
water line and the ks2el, with width increasing in the 
vertical direction according to an average deadrise 


angle. 
The average forward sidewall area is found to be: 
(see fig. 2 and 3b) 


AVERAGE WIDTH = LDBA2R1/(2*TAN(DR1)) +*4S10 


СЕТЕ 


whare LDBAR1 is draft at forward centroid, DRI = 
average deadrise angle forward and 4510 15 th2 average 
Keel width of tne forward section. Therefore submerged 


volume of a forward sidewall is given by: 


VOL = L1*(LDBAR1/(2* TAN (DR1)) +¥S10) *LDBAR1 


(ТТ1-2) 


where L1 is the length of a forward buoyant section. 


Following the same procedures tha after sidawall width 


16 





"s 


(see Dod апа 92) 


AVERAGE WIDTH = LDBAR2/(2*TAN (DR2)) +WS20 


(II*3) 
where DR2 = average deadrise angle aft 
LDBAR2 = draft at after buoyant centroid 
WS20 = average keel width aft 
Submerged volume of an after sidewall 15 then 
determined to be; 
VOL = L2* (4S20+LDBAR2/(2*TAN(DR2))) *LDBAR2 
(II-4) 


where L2 is the length of the after sidewall. 


т 
№ 


. SEAL MODELING 


A major change іп McIntyre's modal was made by 
remodeling the stern 5231. In tha original model a large 
negative moment was required to place the pitch 
equations in steady state equibibrium. Investigation of 
seal forces as presented іп (91 revealed a large 
discrepancy between the after seal forces generated іп 
McIntyre!s nodel and calculations based оп the 
information in (81 . This discrepancy was further 
supported Бу actual s3al load data presented by Layton 
in Ud . All seal dimeasions used in this nodel were 
taken from data prasented in ПП . Seal hinge Locations 


were taken from the data used in GI. 
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Seal lift force forward is generated by plenum gauge 
pressure acting downward on the seal face area in 
contact with the water passing underneath the seal. 
This downward force results in a lift force being 
generated as the downward movement of the seal iisplaces 


the water beneath it. 


The seal is considered to have a straight forward face 
Which runs from the hinge point at the plenum roof to 
the bottom of the k3el at an angle determined to be 
best-fit from the data in ref. OU . The wetted surface 
of the seal is considered to lay parallel хо and on the 
water at the point where the seal face intercepts the 


calm water line. as a result forward seal lift equations 


are: 
ASEALI = WIDTH*XSEALI 

(II-5) 
XSEAL1 = (LD-L3*TAN (THETA) ) /SIN(31) 

(II-6) 
LD = draft at the C.G 
Б - Gisctance from Z.G to point of contact for forward 


seal race. 
LD ~ L3*TAN(THETA) = draft at forward seal face contact 
point 


31 degrees is the interior angle of the seal. 


The after seal wettel area 15 determined in a similar 


manner to pe 


ASEAL2 = WIDTH*XSEAL2 


18 





(DT 


XSEAL2 = (LD+L4*TAN(THETA))/SIN(32) 


(TES) 


In this model forward seal pressure is considered 
equal to plenum pressure. As a result forward seal lift 


ls : 


HSSALP = PBBAR*ASEAL1 


(II-9) 


Since plenum guage pressure tends to force the after 
seal upward, resulting in venting of the plenum air 
mass, the rear seal is operated at à pressure slightly 
above that of the plenun to reduce the venting tendency 
of the after seal. As a result the after seal r2guires a 
modification in 00331119. Іп McIntyre's work the 
pressure differential (PSEAL-PBBAR) only was considered 
in the Lift force. However, as previously statad, there 
were large discrepancias in computed loads as compared 
to measured data from [10]. Therefore in this model, 
əd to be PBBAR plus 


differential pressure. AS a result after seal lift is 


er 
the seal working pressure is consider 


given by: 


HSEALA = (PBBAR+PDIFF) *ASEAL2 


(II-10) 
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Е. PLENUM MODELING 


In the work presented by McIntyre the plenum was 
considered a box with the only dynamic dinensional 
Variation being in the Z direction. However (as can be 
seen in fig. 1), althaagh the forward edg2 can be 
assumed to be reasonably constant in tne Y-Z plane, the 
forward face of the rear seal slopes rather steeply. 
Therefore the plenum length (in the X direction) varies 
rapidly as draft varies. As a result the 2ffective 
plenum roof area consists of that area from the after 
edge of the forward seal to the forward hinge point of 
the rear seal, plus the vertical projection of the 
unwetted forward face of the rear seal. AS an 
additional consequence it can be seen that the effective 
plenum center of pressurz= varies fore and aft as draft 
varies . In fact as draft increases the C.P. moves 


тогчага. 


The following variable names will be used to describe 


the plenum chanber (ses fig. 4.) 


LPR = length of plenun roof from the rear 2132 of the 
forward seal to the hing? point of the after seal. 

LPWL = Length of plenum at the keel from the after edge 
of the forward seals to the transom. 

LP = length of plenum at the waterline. 

BUBHGT = vertical distance from the plenum roof to the 
bottom of the keel at the center of gravity. 


AB = plenum effective roof area. 


Using the previous d2finitions, the plenum roof length 


and area is: 
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LP = LPWL~ (LPWL-LPR) *(LO/BUBHGT) 
(II- 11a) 


and effective plenum roof area is 


AB = WIDTH*LP 
(II-11b) 


with the result that lift due to plenum gaug2 pressure 


Ec. 


HPRES = ~AB*PBBAR 


СЕЕ hic} 


As a comment оп the complexity of the Lift forces of 
this model note that for constant plenum pressure 


increased draft results in decreased plenum lift. 


G. CRAFT DYNAMICS 


1. Plenum air Mass ani Pressure Equations 


(= | ——— > => -h p e 


For the purposes of this model the rear seal 
leakage area is considered constant since to assume 
Otherwise would lead to modeling a very difficult and 
poorly understood dynamic process. Іп addition the 
plenum pressure VS air nass process is assumed to be 
adiabatic in nature. Pan characteristics are those used 
in the constant air leakage heave only model >f Gerba 
and Thaler іп (ГІШ. Air flow into the plenum is taken 


as positive. 
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Therefore the aicflow into the chamber is given 


by: 


Qin = N*(QI0-Kg*PBBAR) 


CLE По) 


Where N is the nuaber of fans supplying air to the 
plenum, QIO is the equilibrium air flow rate for the fan 
With zero gauge pressure and Kj is the slope of the 


fans! airflow VS gauge pressure curve . 


Qout is airflow out of the plenun and 


н 
и 


determined py the seal leakage equation: 


Зоц= = Cn*Al*(2*PBBAR/RHOA)**1/2 


(II-13) 
Where Сп = seal Leakage coefficient, Al = seal leakage 
area, and RHOA = mass density of air. Therefore nass 
flow rate is: 
“BDOT = RHOA*(Qin-Jout) 
(II-14) 
The plenum air dynamics are considered in 


equilibrium when Qin = Jout 


py 


Absolute Plenum Pressure is etermined from the 


adiabatic gas law to be: 


Pb - Pa*(MB/(RHOA*PBBAB) ) **GAMMA 
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(ITAL) 


and, 


PBBAR = Pb-Pa 
(II-15a) 


Where Pa = ambient atmospheric pressure (LBF/FI**2 GAMMA 


= adiabatic process coefficient, Pb = plenum absolute 
pressure, and PBBAR = plenum Jauge pressure. 
2. Buoyant Forces 





Buoyant force is calculated 27 the submerged 
volume multiplied by the mass density of the fluid times 


ENS gravitational constant. 


Therefore forward buoyant force is: 


dB? = -2*RHO*G*VOL 


(II-16) 


Where volume is given by eqn. II-2, Rid is the 
mass density of water (slugs/ft**3) and з is the 
gravitational constant. ЕО ОЕ 2:15 introduced 
to account for 2 sidewalls. The total expression for 


HBF is: 


HBP = -2*RHO*G*L1*LDBAR1* (LDBAR1/ (2*TAN (DR1)) +WS01) 


(11-17) 


Following the same pattern, 
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НВА = -2*RHO*G*L 2*LDBAR2* (LDBAR2/(2*TAN (DR2)) =W510) 


(11-18) 


Forward seal force is determined from eqn. II-9 to 
be: 


HSF = ~PBBAR*WIDTH* (LD~L3*TAN (THETA) ) /SIN (31) 
y 


(1т- 19) 


where HSF is the heave force due to forward seal. 


In a similar manner HSA is: 


HSA = (PBBAR+PDIFF) *WIDIH* (LD+L4*TAN (THETA) /SIN (32) 


(ІІ-20) 


McIntyr2 introduced the planing force phenomena in 
his model to account for the proven ract that the CAB 
SES tends to decrease it's draft and pitch down as speed 
increases. His development of tha planing equation was 
not well explained aad the basis for selection of the 
planing coefficient was not clearly stated. Raidel [12] 
did a special investigation or the planing force 
phenomena. His results were based on flat plate planing 
effects to approximate the flat bottom of the craft 


Sidewalls. 


24 





Since the water is in contact with the vertical 
portions of the sidewalls it was felt a more accurate 
representation was neeced. Barnaby (0131 states that 
water contact on the vertical sides of a planing hull 
rapidly reduces planing lift advantage since one of the 
advantages of a planing hull is reduced drag due to 
wetted surface. planing lift results in a speed 
sensitive tradeoff between static lift (buoyant effects) 
and dynamic lift (planing effect) that іс directly 
proportional to the angle of attack (for smali planing 
angles) and the square of the velocity.  Theretore the 


planing force is of the form: 
HPLAN = K*V**2*THETA 


From Reidel's work, 


HPLAN = 0.5*AREA*V**2*PI*SIN(THETA) 
(II-2 1) 


for small angle approximations tnis siuplifies го: 


HPLAN - RHO*AREA*V** 2*PI*THETA 


Which is the planing force for a flat plate. Barnaby 
gives no simple formula for the losses due to inmersion 
so we are left with an arbitrary loss coefficient, 
PLCOEF such that: 


HPLAN - (PLCOEF) *2X*RHO*AREA*V**2*PI*THETA 
(II-22) 


Phe added factor of 2 accounts for 2 sidewalls. 


The planing area is taken from the scale drawings 
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developed from sidewall data, and the force is assumed 
to act at the centroid of the flat keel surface, which 
is consistant with Reidel's work. The loss coefficient 
will ре developed in the section on achieving 


equilibriun. 
5. Plenum Gauge Pressure Lift 


Combining eguatisns II-11a ,II-11b0b , and 1:-112 


results in: 


HPRES = -PBBAR*WIDTH*(LOWL- (LPWL-LPR) * (LD/BUBHSD) ) 


(11-23) 


Lift Total = HPRES + HB& + HBA * HSF * HSA + HPLAN 


4hich must exactly equal craft weight (W) £or the hull 


во Бе im vertical equilibrium, or: 


W + HPRES + HBF + HBA + HSF + HSA + HPLAN = 0 
(II-24) 


defines 2quilibrium condition in heave motion. 


A RAFT PITCH MOMENT EQUATIONS 
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Forward buoyant moment 1s the forward buoyant lift 


times the effective lever arm, Or; 


PBF = -HBF* (Lever arm forward) 


-ШВЕ ТВО 


РВЕ 


(11-23) 


Note that the negative sign results from the 
definition of positive noments and the fact that upward 


forces are negative by iefinition. 


i 


PBA HBA*LEVER AR 


HBA *L6 


PBA 
(ІІ-26) 


using the same reasoning, after buoyant moment is: 


BE total = PBF + PBA 


(II-27) 


D Seal Moments 


= җе с=ш= сезщ Cn ccm amen oer se cee 


Since the seals arce relatively distant from the 
C.G. dynamic variations in the lever arms are considered 


in seal moment equations and are developed ners. 


27 





The seal тоге is modeled tc act at the 
longitudinal center ofthe seal wetted length. Therefore 


the effective lever arn for the seal moment is: 


ВЕБЕ = L3 — XSEALI/2 


where PLSF - moment lever arm for seal forward. Ас а 


result: 


PSF = -HSF*PLSF 
PSF = -HSF*(L3-(XSEAL1/2)) 
(2) 
In a similar manaerc 
PLSA = L4 + XSEAL2/2 
PSA = ESA*FPLSA 
PSA = HSA* (L4+XSEAL2/2) 
(11-29) 


PS TOTAL = PSF+PSA 


Note that the dynamic pitch sensitivity to moment 
arm length is included in the XSZAL tarms as defined in 


equations II-5 and II-3. 
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For the planing force developed in equation II-22 
the planing force is concentrated at the centroid of the 
flat keel which results in: 


PPLAN = L7*HPLAN 


(II-30) 


First the dynamic variation in effective center of 


pressure will be developed. (see fig.(4)) 


plenum length is taken from equation [I-11a to be 


LP - LPWL- (LPNL-LPR) *(L2/BUBHGT) 


Berining LDIFF to be 


LDIFF = LPWL-LER 


(II-31) 
and, 
XCP = XCPO+ (LDIFF*LD/BUBAGT) /2 

(ІІ-32) 
where 
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ХСРО = LPWL/2-XCG 
(II-33) 


and ХСС is the distance from the transom to the C.G. 


(see fig. 4.) Finally, 


PPRES - -HPRES*XCP 


(II-34) 


5. Pitch Daaping 


Any frictional damping terms dependent upon the 
Square of pitch rate (1.2. tangential velocity about the 
Y axis) disappear when the Taylor Series expansion is 
carried about the zero rate equilibrium point. Since 
for small signal approximation the model should reflect 
actual craft charactecistics, it was hoped that the 
changes in seal and plenum modeling would reveal the 
Missing damping characteristics. As it turnei out this 
was not the case. Examination o£ the hydrodynamic 
equations in [9] revealed an added mass effect dependant 
upon a cross velocity term between linear X direction 
velocity U, pitch rate Q, and added-mass at the stern 
A33S due to slender body theory apolied to the sidewalls 


. This moment is: 


FP = -A33S*XLSS**2*ÜX*THETADOT 
(11:35) 


Where XLSS is the distance from th2 C.G. to the stern, 
and THETADOT is pitch cate. 
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В = SS **2 Y RHO*PI/S 


(II-36) 


where YSS is the cross section width of the side 
Wall at the stern waterline. Therefore, 


FPDAMP = -2X*XLSS**2*A33S5*U*THETADOT 


ZT DKRUTFRL IS RAD STIL SSH 2SRHOFUFTHELDOT ИЗ 
(ез 


Miert actor ЭЕ 2 LS introduced to account for two 


sidewalis. 


Since the moments are all referanced to the craft 
center of gravity, the craft weight contributes zero 
moment in that frame ог reference, Therefore total 


moments are: 


PSF + PSA + PBF + PBA + PRES + PPLAN + PDAMP = 9 


(ІІ-38) 


І. HEAVE ACCELS£RATION EQUATIONS. 


From Newton's law, acceleration and force are related 


DY: 
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FORCE = MASS*ACCELERATION 


For the linearized model we are interested in the 


acceleration on the craft and as a result: 


ACCELERATION = РОВСЕ/ MASS and therecore the heave 


equation of motion becomes: 


ZDDOT = (HBF+HBA+HSA+HPLAN+HPRES) /MASS 
(ET=3) 


where ZDDOT is the second derivative of Z witı Tespect 


to tine. 


J. PITCH ACCELERATION EQUATIONS 


Angular acceleration is: 


THETADDOT = PTOTAL/IYY 


Е ТҮҮ Іс -hə Сеа те inSetial momant about tas Y axis. 


As a result total angular acceleration is; 


THSTADDOT = (PBE+PBArPST+PSA+rPPLANFPPRES+PDAHP) /LIY 


(II-40) 


This completes the derivation of all aonlinear 
=quations needed to describe the two-degree-o0f£-freedon 
MAD SES. In the following chapter these aodnlinear 
2quations will be linearized using the Taylor Series 
expansion about a steady state operating poiat which 


will result in the equations needed to generate a 
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dynamic digital-computsc simulation of craft notion in 
the time-domain. 
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LII. LINEARIZATION OF THE ASAVE-PITCH EQUATIONS. 


A. TAYLOR SERIES EXPANSION 


For small perturbations about a steadily state operating 


point the nonlinear function Z can be approximated by: 


EZ E, Y) 


2 {0) +92 = Z(X(0) ,Y(0))+ (92/9X) *d42+ (92/9Y) *dY 


where (23Z2/8X) and (32/9Y) are the partial derivatives 
of the function 2 evaluated at the steady state 
Operating point X(0),Y(0). Cancelling the stealy state 
values from both sides of the equation yields the 


differential model equations: 


47 - (97/ӘХ) х4х%(97/2Ү) «ау 


ада 


2 = 2(0)+47 


where dz,dX, and dY are tne dirferential variables for 
Ze X, Y respectively, (02/9Х) апа (dZ/dY) are the 
partial derivatives of Z2 evaluated at the steady-state 


operating point. 
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В. LINEARIZATION OF THE MODEL EQUATIONS 


There are basically four equations to linearize for 


the two-degree-of-freedon model: 


ZDDOT = F(Z,THETA, MB) 


THZTADDOT = G(Z, THETA, THEPADOT,MB) 


MBDOT Н (2 „МВ) 


| 


P3BAR - І(2,МВ) 


therefore the equations for the differential state 


variables are: 


1ZDDOT = (4F/d4Z) *dZ+ (dP/ATHETA) *ATHETA+ (AF/dMB) *амв 


(III-1) 


ITHETADOT = (46/427) «47 + (16 /4ТНЕГА) «ЗТНЕТА + (Ус/ МВ) *d MB 


(218-2) 


11300т < (4Н/42) *4Z2 +(11/4THETA) *dTAETA+ (dH/AMB) *AMB 


(III-3) 
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APBBAR + (471/42) «47 + (4Г/АТНЕГА) ЖаТНЕТА + (ат /дМв) «див 


(III-4) 
Note that the left ha nd side variables are 
differentials, not time derivatives. 
Ва LINEARIZED EQUATIONS 
1. Plenum àir Mass Flow Derivatives 
Qin - N*(QiO-Ki*PBBAA) 
(GIST) 
Qout = Cn*Al* (2*PBBAR/RHIA) **17/2 
(11-13) 
MBDOT = RHOA*(Qin-Qout) 
(II-14) 


as a result the differential equation for MBDOI becomes: 


амарот = RHOA*(dQin-dQout) 


(III-5) 


dQin = (dQin/dZ) *dZ+ (dQin/dPBBAR) *dPBBAR+ndQin/i4B) жамв 
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(ITI-B) 


dQin/dPBBAR = -N*K3 


The partial derivatives of Qin with respect to all other 


variables are zero. 


dQout = 
(dQout/dZ) *dZ+dQout/dPBBAR) *dPBBAR+(dQout/dMB) *dMB 


dQout/dPBBAR = Cn*A1L*((2RH0A/2* PBBAR) **1/2) /RHOA 
(ІІІ-7) 


The partial derivatives of Ооп] with respect to all 
other variables are zer. Inserting equations III-6,and 
III-7 into 11-14 у1е1 45 the eguation СОР the 
differential variable dMBDOT: 


dMBDOT = (-RHOA*N*KI-Ca*AL*((RHOA/2*PBBAR) **1/2)) ХаРВВАВ 


(11-15) 


2. plenum Gauge Pressure Derivatives. 


Pb = Da*(Mb/RHOA*VD) «“ЗАММА 


(ТЕЛ) 


ana 


PBBAR = Pb-Pa 


S 





аса Ра 15 considered constant within the time frame of 


this model, 


dPBBAR = dPb 
(ILI-9) 
As a result, 
dPBBAR = 
GAMMA* Pax (MB (0) /RHOA*Vb (0) ) ** (GAMMA- 1) *d (UB/RHOA*VB) 


where 


4 (MB/RHOA*VB) = (Т/ВНОАЖУВ (0) ) *4MB+ (MB/RHOA) *4(1/VB) 


d (MB/RHOA*Vb) (1/RHOA*VB(0)) *4MB-(MB/BHOA*V3**2) *dVB 


d (VN-AB*LD) 


И 


ачв 


dave = -АВ(0)хаһ,о/42-10(92) *dAB/dZ 


remembering that: 


110/42 = 1 


we now need the partial of efrective plenum roof 


area with respect to Z. Defining 
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AB = ABWL-(ABWL-ABR) *LD/BUBHGT 


(ІІІ-10) 


ABWL = WIDTH*LPWL 


ABR = WIDTH*LPR 


И ече derining 


ADDIFF = ABWL-ABR 


which іс а constant, tha effective plenum roof area is: 


AB = ABWL-ABDIFF*LD/BUBHAST 


using the previous definitions the partial of 2ffective 
roof area becomes: 
dAB/dZ = -ABDIFF/BUBHGI 

(ІІІ-11) 
and the differential volune 15: 
dVb =-(AB(0)-LD (0) *2*(ABDIFF/BUBHGTD) ) *dZ 

(ІІІ-12) 


Rich the result that 


4 (4B/RHOA*Vb) = (1/RHOA*VD(0)) *dZ + 


(Мо (0) /RHOA*Vb (0) **2) * (AB (0) -2*LD (0) *(ABDIPF/BUSHGT) ) *dz 
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Finally, 


dAPBBAR = GAMMA*Pb(O) *((1/Mb (0) ) *iND + 


(1/Vb(0) * (AB (0) -2*LD (0) * (ABDIFF/3UBAGT) )) *4Z 
(ІІІ-13) 


The total heave force equation is: 


H(total) = HBFtHBAtHSFtHSAt+HPLAN+HPRES 


The total partial derivatives will be done term by 


tern. 


For the forward buoyant force the following 
simplifications are made: 


HBF » -2*BRHO*G*LÍ*LDBART* (LDBAR1/(2* TAN (DR1)) *W310) 


let К1 = 2*RHO*G*L1 then, 


HBP = -K1*LDBAR1* (LDBAR1/ (2*TAN (DR1}) ) +¥S10) 


where LDBAR1 = (LO-LSTAN (THETA)). Combining the 
previous two equations, 


HBF = ~K1*(LD-15*TAN (THETA) ) * 
((LD-L5*TAN (2 HETA)/ (2*TÀN(DR1) ) * S10) 


40 





assuming small Deen angles, TAN (THETA) can be 
approximated by THETA which results in the following 


Simplification: 


HBF = 
-K1*(LD-LS*THETA)*((LD-L5*THETA) / (2*TAN (DR 1) ) *i310) 


Taking the partials of ABF with respect t> 7 апа 
THETA yields: 


d(HBF) /dZ * -K1* ((LD (0) -L5*THETA (0) /TAN (DR1) *3510 
(177-14) 


апа 


d (HBF) /dTHETA = 
K1*L5* ((LD(0) -LS5*THETA (0))/(2*TAN(D/r1))+WS10) 
+K1* (-L5/(2*TAN (DR1)) * (LD(0) -LS5*THETA (0)) 


Which Simplifies to : 
d (HBF) /4ТНЕТА = 
K1*L5* ((LD (0) -LS*THETA(0)) /(TAN (DR1) ) +9510) 


(ІІІ-15) 


Using the same approach, the partials for buoyancy 


force aft are: 


4 (HBA) /dZ - -K2* ((LD (0) *L6*THETA (0) ) /(TAN (DR2)) * $S20) 


(III-16) 
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d (HBF) /dTHETA = 
-L6*K2* ((LD (0) +L6*THETA (0) / (TAN (DR2) ) +WS10) 


(III-17) 
For the forward seal the total differentíals are: 


DHSF =(d(HSF) /dZ) 42 + (4 (Н5Р) /АТНЕТА) жаТНЕТ 
+ (9 (Н5Е) /АРВВАВ) *АРВВАЕ 


the following simplification is made; 


HSF = -PBBAR*WIDTH* ({LD-L3*TAN (THETA) /SIN (31) 
(II-19) 
let WIDTH/SIN(31) = K3 and TAN(THETA) = THETA, then, 


вне —KI~PBBAR* (LD=L3* THETA) 


with the resultant partials being: 


d(HSF) /dZ = -K3*PBBAR(0) 

(III-18) 
1(HSP) /ATHETA = K3*L3*PBBAR(O) 

(ІІІ-19) 
3 (Н5Е) /АРВВАВ = -КЗ* (10 (0) -L3*THETÀ (0) ) 

(ІІІ-20) 
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The partials of after seal force are founi in the 


Same manner to be: 


ан5А = (4 (HSA) /dZ) *GZ+(d(HSA) /dTHETA) *4 THETA 
+(d(HSA/dPBBAR) *dPBBAR 


HSA = - (PBBAR+PDIPF) *WIDTH* (LD+L4*TAN (THETA) /SIN (32) 
(II-20) 


let K4 = WIDTH/SIN(32) and assume small THETA so that; 


HSA = -KU* (PBBARtPDIFF) * (LD+L4*THETA) 


the partials then become: 


d (HSA) /dZ = -K&* (PBBAR()) +PDIFF) 

(ІІІ-21) 
d (HSA) /ATHETA = -K4x* (PBAR (0) +PDIFF) «14 

(ІІІ-22) 
А (ЕЗА) /АРВЗАВ = -Кц* (1.0 (0) +14*ТНЕТА (0) ) 

(III-23) 


Next tne planing force partials are derived: 


МЕБАП = —-2*XPLCOEF*RHO*AG3EA*V**2*XDpDI*THETA 


let KPLAN = (2*2LCOEPF*RHO*AREA*PI) then, 
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HPLAN = -KPLAN *V**2*THETA 


and, 


а (НРІАМ) = (d(HPLAN) /dTHETA)*dTHETA 


for this constant speed model, therefore, 


3 (НРЬАМ) /АТНЕТА = -KPLAN*V (0) **2 
(III-24) 
Plenum Gauge Pressure lift was defined to be: 
HPRES = -AB*PBBAR 
(II-11c) 
as a result: 
і (НРВЕЅ) /аРВВАК = -АВ (0) 
(ETT) 
4. Pitch Moment Partials 
The total moments acting on the craft are: 
PTOTAL = PBF+PBA+PSFHPSA+ PPLAN+PPREStPDAMP 
(II-38) 


The partials of II-38 will be found tern by tern, 


peginning with the forward buoyant moment: 
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Par - —HBF*L5 


(II-23) 
therefore: 
d(PBF) = -L5*d (HBF) 
so that: 
d (PBF) = -L5*(d (HBF) /d2) *dZ-L5*( (4d (HBF) /Я4ТНЕТА) «ЯТНЕТА 
(ІІІ-26) 
and the partial of the after buoyant moment is: 
PBA = LO*HBA 
(ІІ-26) 
аРЗА = Го жанВА 
tharefore, 
4 (РВА) = 16* (а (НВА) /42) *12+16* (4 (НВА) /АТНЕТА) «АгНЕТА 
(ІІІ-27) 


Seal forces partials must include thə dynamic 
variation in lever arm length, which for the forward 


seal 15: 


Dom - —HSFP*PLSF 


D. 
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PSF = (PBBAR *ASEAL1) *PLSF 


therefore, 


d (PSF) /GPBBAR = ASEAL1 (0) *PLSF (0) 
(ІІІ-28) 


апа 


d (PSF) dZ = PBBAR (0) *d(2LSF*ASEAL1) /dZ 


where 


á(ASEALT*PLSF)/dZ s PLSF(0)*d(ASEAL1)/àZ + 
ASEAL1 (0) *d (PLSF)/1Z 


ASEALI = WIDTH* (LD-L3*TIETA) /SIN (31) 


and 


PLSF = L3-(LD-L3*THETA) /2*SIN (31) 


therefore, 


4 (2152) /42 = -1/(2*SIN(31) 


апа 


d(ASZAL1) /dZ = WIDTH/SIN (31) 


as a result, 


d (PSP) /dZ = PBBAR(0) *((PLSP (0) *WIDTH) /SIN (31) - 
ASEAL1 (0) /(2*SIN (31) ) 
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(ттт=29) 


а (рѕЁ) /ATHETA = PBAR (0) *d (PLSP*ASEAL1) /dTHETA 


d (PLSF*ASEAL1) /dTHETA = 
PLSF (0) *d (ASEAL1) /ATHETA+ASZAL1 (0) *d (PLSF) /ATHETA 


d (ASEAL1)/dTHETA = -L3*AHIDTH/SIN(31) 


а (РЪЗР) /АТНЕТА = L3/(2*SIN (31) 


as a result, 
1 (BSF) /ATHETA = -PB83AR (0) * ((PLSF (9) *L3* WIDTH) /3IN (31)- 
ASEAL1 (0) *L3/ (2*SIN (31) ) 


(III-30) 


For the after seal the derivatives procesi in the 


manner: 
PSA = -(PBBAR + PDIFF)*ASZAL2*PLSA 
d (PSA) /APBBAR = -ASEAL2 (0) *PLSA (0) 
(LIL > 31) 
In the same manner as shown іп the previous 


paragraph, the partials of moment due to seals aft with 


respect to the other variables are: 
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d (psa) /dZ = -РВВАВ (0) “ ( (РОЗА (0) *WIDTH) /SIN(32) + 
ASEAL2 (0) / (2*SIN (32)) 


(EIT-32) 
d(PSA)/dTHETA = -PBAR(O)*((PLSA(O0) *LA*WIDTH)/SIN)32)* 
ASEAL2(0) *L4/(2*SIN (32)) 
(III-33) 
The partial of planing moment is: 
PPLAN = L7*HPLAN 
(11-30) 
d(PPLAN) = L7*d(HPLAN) 
4 (PPLAN)/dTHETA = -L7*KPLAN*V (0) **2 
(LII-34) 


For this constant speed model all other partials are 
identically zero. The Plenum Gauge Pressure also has 
dynamic variation in the moment lever arm and the 


derivatives are: 


> —ХСР*НРКЕ5 


(11-34) 


therefore, 


APPRES = -НРВЕЗ (0)4 (ХСР) -ХСР (0) «4 (НРВЕЗ) 
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where 


XCP = XCPO+ (LDIFF*LD)/(2*BUBHGT) 


therefore, 


d(XCP)/àZ = LDIPF/(2*BUBHGT) 


апа 
а(РРВЕ5) = -((HPRES (0) *LDIFF) /(2*BUBHGT) ) *dz 
-XCP(0) *d(HPRES) 


substituting equation (III-25) yields: 


4 (РРВЕЗ) = -((HPRES(0) *LDIFF) /(2*BUBHGD) ) *dz 
*XCP (0) *AB (0) *dPBBAR 


Finally, the partial of pitch dauping moment is: 


PDAMP 2 -2*PI*XLSS**2*YLSS**2*nRHO*u*THETADOT/S8 
(II-37) 


let 


A2*PI*XLSS**2*YSS**2x*xRHO/8 = KDAMP 


then 


PDAP - -KDAMP*V*THETADOT 


aS a result, since this is a constant speed model, 
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3 (РОАМР) /АТНЕТОГ = -КОАМР*У (0) 
(ТТТ=35) 


where THETDT is the first derivative of theta with 


respect *o time. 


5. Total Differentials 


Since F = M * A and the interest ls in 
d(F)/M. the result 
being that all differentials пузо ое divided by craft 


H 


acceleration terms, A = F/M and dA 


mass to achieve the jdesired linearized heave results. 
The same general comment holds true for moment 
differentials, that is; they all must be divii2d by the 
craft moment of inertia about the Y-axis. 


All differentials must be sorted by the variable 
of differentiation so that the Linearized differential 


equations of motion may be written. As an example; 
AZDDOT = ((d (HBF) /dZ) + (4 (ЧВА) /42)) /МА55) х42 + 


The State Variable representation of tha system 
Will be presented in a later section of this chapter, 
however it is convenient to define the State Variables 


at this time. 
8 
7, БОТ, СТНЕТА, ЗРНЕТОГ, Я МВ) 


(ІІІ-36) 


апа 
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У = iPBBAR 
(ІІІ-37) 


since PBBAR is not a state variable the definition for 
dPBBAR must be substituted for dPBBAR in all equations. 
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Here the total differentials will be found. No 
2xplanatory taxt Will Бе included unless special 


substitutions requiring an aside are necessary. 


The buoyancy term3 are: 


DHBZ = d(HBF)/dZ + 4 (НВа) иат, 


(грива) 


A term must be included which was not previously 
derived to account for the sensitivity of seal forces 


даз to PBBAR since it is not a State Variable. 


DESPZ = (4 (ESF) /APBBAR) * (d (PBBAR) /dZ) 


as a result, 


4 (Н5Е) /4РВВАЕ = -K3* (LD(0) -L3*THETA(O)) 

(ПІЕЕЕ- 20) 
and 
DHSFPZ = -K3*(LD(0) -L3*THETA(O)) *d (PBBAR)/dZ 


sinilarly foc the after seal moment: 
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d (HSA) /APBBAR =-K4*(LD(9) +LU*THETA (0)) 


(III-23) 
and 
DHSAPZ + - Кж (10 (0) *LU*THETA(0)) *d (PBBAR) /dZ 
.The total seal force differential in Z becomes: 
DHSZ = d(HSF)/dZ + dftHSA) /dZ +DHSFPZ + DHSAPZ 
(ІІІ-39) 


For plenum gauge pressure lift differential the 


Same substitution must be made. 


ddPRES = -AB(0) *dPBBAR 


therefore, 


DHPZ = -AB(0) *d (PBBAR) /4Z 


(ІІІ-40) 


Define the sensitivity coefficient of 2000Т Хо 7 
to be: 


DZZ = (DHBZ + DHSZ + DHPZ) /MASS 


(III-41) 
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The following partials are now defined: 


ОНЗТН = 4(НВЕ)/АТНЕТА + d(HBA) /ATHETA 

(III-42) 
DHSTH = d(HSF) /ATHETA +1 (H3A) /ATHET 

(III-43) 
DEPTH = d(HPLAN) /ATHETA 

(III -44) 


Define the sensitivity coefricient of ZDD3T due to 
УНЕТА Со Des 


DZTIH = (DHBTH * DHSTH * DHPIH)/MASS 


(РЕ) 


4. Total Differentials due to 43 


dere since dPBBA3 is not a state variable the 
partial of PBBAR with respect to MB must be substituted 


in the plenum lift and seals lift terms. 


à (HSF) /dPBBAR - -K3* (LD (0) -L3*THETA (0) ) 
(III-20) 


Dy substituting d(PBBAR) /dMB for dPBBAR the result is: 


S 





d (HSFP) /AMB = -КЗж (10 (0) -L3*PHETA (0) ) *à (PBBAR) /1M8 
(III-46) 


similarly, 


-Куж (10 (0) + ЦТ НЕТА (0) ) «4 (РВВАВ) /44B 


d (HSAP) /dMB 


(III-47) 
d (PRES) = -AB(0) *dPBBAR 
using the previous substitution yields: 
d (APRES)/d4B = -AB(0)*1 (PBBAR) /dMB 
(IIT-948) 


The sensitivity coefficient of haave due to MB can 


now be defined to be: 


0248 = (d (HSFP)/dMB *d(3SAP)/dMB * i(HPRES)/dMB) /MASS 


(III-49) 
3. Total Pitch Moment Differentials ia 2 
DPBZ = d(PBF)/dZ + 4 (РВА) /47 
(III -50) 


The dPBBAR substitution yields the partial of 


seais with respect to pressure with respect to Z. 


54 





ASEAL1 (0) *PLSF (0) *d (PBBAR) /dZ 


DPSFPZ 


(III-51) 


DPSAPZ = -ASEAL2(0) *PLSA (0) *d (PBBAR) /dZ 
(ІІІ-52) 


therefore the total seal differential in Z becones: 


DPSZ = DPSFPZ + DPSAPZ “4 (Р5Ғ)/42 ж 4 (Р5А)/42 


(III-453) 
Next the Plenum 3aige pressur2 tarü is isrived: 


4(PPRES)/dZ = -HPRES (0) *LDIFF/(2*BUBHGT) 


4 (РРАЕЗ) /АРВВАЕ = АВ (0) *ХСР (0) 


therefore, 


DPPRESZ = АВ (0) *ХСРр (0) «а (РБВАБ) /27 


and 


DPRESZ = d(PPRES) /dZ +d(PPRESZ) /iZ 


(ІІІ-54) 


The sensitivity coefficient of THETA with respect 


to Z becomes: 
DTHZ = (DP5Z + DPSZ * DPPRES)/IIY 
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(III-55) 


10. Total Moment Differentials in IHETA 
DPBTH = (d(PBF)/dTHETA + d(PBA) /ATHETA) 
(ІІІ-56) 
DPSTH = (d(PSF)/dTHETAÀA * d (PSA)/dTHETAÀ 
(III-57) 
DPPLTH = d(PPLAN)/dT HETA 
(ІІІ-58) 


Cri cleanc. core Panta tobe. 


а 
O 
(0 


define the sensitivity 


DTHTH = (DPBTH + DPSTH + DPPLTH) /LYY 


(TELS) 


Also, define the sensitivity coetfficient in THETDT 
to be: 


DTHDTH = (d (PDAMP) /ATHDT) /IYY 


(LII-60) 





dPSAMB = (d(PSA) /APBBAR) * (d (PBBAR) /dMB) 

(III-61) 
АР5ЕМВ -(4(Р5Ғ) /4РВВАВ) % (d (PBBAR) /AMB) 

(ІІІ-62) 
DPSMB = DPSAMB + DPSFMB 

(III-6 3) 


Inserting the dPBBAR substitution yields the partial of 


plenum gauge pressure moment with respect to MB: 


а (РРНЕ5) /АРВВАВ = -ХСР (0) «DHPRES 


4 (НРВЕЗ) /АМВ = -АВ (0) *3 (РВВАВ) /dMB 


ОРРМВ = ХСР (0) “АВ (0) *d (2BBAR) /dMB 


(III-64) 


Define the sensitivity coefficient of THETA due to 
MB to be: 


ОТНУВ = (DPSMB + DPPM3)/IYY 


(III-65) 
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4 (4В) /4МВ = (d(MB) /dPBBAR)* (d (PBBAR) /dMB) 


define the sensitivity coefficient to be: 


DMBMB = 4 (МВ) /дУВ 
(ІІІ-66) 


define the sensitivity coefficient in Z to be: 


DM3Z = d(MB) /dz 


(ІІІ-67) 


DE ZETOTAL HODEL SENSITIVITY COSLSPPILCIENTS 


1. Heave Sensitivity Coefficients 


The sensitivity coefficients for the linzarized 


dizferential equation in heave acceleration are: 


DAZ, 20400, and 02 МВ. 


2. Pitch Sensitivity: Coeffi 


a. 


The sensitivity coefficients for the linearized 


differential equation in pitch acceleration are: 


DTHZ, DTHIH, DTADTH, and DTHMB. 
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The sensitivity coefficients for plenum air mass 


are; 


DMBZ, and DMBMB. 
4. Plenum Gauge Pressure Sensitivity Coefficients 


. The sensitivity coefficiants for Planun Gauge 
Pressure are: 
DPBZ, and DPBMB. 


E. STATE SPACE REPRESENTATION 


The state matrix form of a dynamic system is: 


XDOT = АХ + BU 
where X = a vector of variables representing the states 


of a set of first order Linear differential 2yuations, 
and U is the vector of forcing functions, and A and В 
are the matrices of the respective gains. The form of 


the equations is: 


X1DOT - X2 
X2D0T - X3 
ESDODOSORI*XT + KZEX2 + K3*X3 
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for the linearized model of the CAB SES the state 


variables were defined to be: 


+ 
X - (42, 42р0Т, THETA, dTHETDT, dB) 


Ав а result the state matrix representation of the 


linearized unforced Heave-Pitch nodel is: 





X 1DOT 0 1 0 0 0 

X2DOT 022 0 DZTH 0 DZMB 

X3DOT | - 0 0 0 1 0 

er DTHZ 0 DTHTH ртнртн ртнив 

X5DOT | 0482 0 0 0 D4348 | 
(III-68) 

and 


dPBBAR = DPBZ*X1 + DPBMB*X5 


#11-67 
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IV. ГІМЕ DOMAIN VALIDATION 


Se AO A ғ жақ u Su сит E E A A 


The linearized model developed in chapter III was 
vaiidated against the t he Six-degree-or-freedon 
nonlinear model. There were several changes made in the 
simulation techniques used, and at the same time the 
changes in center of pressure dynamic variation and 
planing forces were entered into the 6 DIF modal. 
Estabiishment of equibibrium for the linear model was 
accomplished Ot тве using a TI-59 programmable 
hand-held calculator with printer. 


The first change in simulation technique is the 
operating point at which the linear sensitivity 
coefficients were evaluated. In 4cIntyre's model these 
coefficients were evaluated at the initial operating 
point. In this work the coefficients are evaluated at 
the final steady-state equilibrium point since this 15 


customary in Linearization work. 


The second major change in simulation was the method 
used to determine the planing force 1055 co2fficient 
(PLCOEFF) and the establishment of the acbitrary 
location of the initial plenum center of pressure. Іп 
this model all heave fòrces ara well defined (as nodeled 
here) With the exception of planing force lift. Asa 
result, craft weight ani all heave forces 2xcapt that 
due to planing lift ar2 summed and the residual assumed 
to be due to the planing phemonema. This residual is 
then used to determine (off-line) the value of PLCOEFF. 
Note that this method also accounts for any Lift due to 


aerodynamic forces оп the бон that wer2 disregarded. 
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Next all moments are suaned about the С.С. except the 
moment due to plenum pressure lift. The residual moment 
is assumed to be due to plenum pressure. This residual 
is used (again off-line) to establish the initial center 
of pressure. This accounts for any moments not 
previously accounted for (such as aerodynamic forces on 
the bow.) McIntyre used the planing lift coefficient to 
force the model into agreement with the 6 DOF nodel for 
natural pitch frequency and used a large, arbitrary, and 
static moment to force the models into steadv state 
pitch angle agreement. It was falt that using tne 
present method ЗЕ determination would laave tha 
equations in a form that would Бе more amenable to 
extraction of the factors most directly affecting craft 


characteristics as nodeled herein. 


A. LINEAR PROGRAM DESCRIPTION 


User's instructions have been included as apoendix A. 
The program flow can be isolated into four Main 


sections. 


1. First the values for craft weight, draft, pitch 
angle, plenum pressure and speed are entered. TO 
initialize plenum air mass, draft is used to solva for 
initial plenum air volume, which is then used with 
initial plenum gauge pressure to solve equation ІТ-15 
for plenum air mass. Equation [I-12 is then solved for 
Qin. Qout is set equal to Qin and equation [I-13 is 
soived for sea: leakage area. Ail other heave and pitch 
nonlinear equations are solved and their values are 
listed along with residuals of heave aad pitch 
summations. These resiiuals are added to the initial 


values of C.G. and angular acceleration and tne initial 
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values are saved for the time domain simulation section. 


2. The same operatioas are carried out For tae steady 
state values of weight, draft, pitch angle, ani plenum 
pressure. The residuals, forces and moments are again 
listed. 

3. These steady state values are used to evaluate 
the sensitivity coefficients lavelopad in chapter III 


which are also Listed. 


4. Finally the stat? variable solution is sinulated 
using a four-point fixed step Runge-Kutta algorithm and 
listed for the print interval desired. The option is 
included to output either printer plot grapnical display 
Ming che NDS printer plot Subroutines PLOTP or the X-Y 
Versatec plot subroutine DRAWP. The latter should be 
used for final runs only, with PLOTP used on 
intermediate runs. Provisions are made for up to 300 
discrete output points for the time domain tabular 


listing. 


В. SIX-DEGREE-OF-PREEDOM HODSL SOLUTION 


Two changes were made in the 6 DOF model to account in 
а more analytic manner for the observed tandancy of the 


XR-3 to pitch down as speed increased. 


The first change was to remove an 2upirical equation 
used to correct the center of pressur2 (51 and replace 
Dt with the dynamic center of pressura variation 
developed in chapter ГГ. 


The second change «as to include in the sidewall 


Subroutine the nonlinear equation for planing forces, 


«ich the loss coefficient adjusted to mak2 up (аз 
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closely as possible) for lift loss due to the revised 


modeling of the plenum chamber. 


Finally, the 6 DOF model was ran at both 0722 lbs and 
6050 lbs weight to establish the initial and final 
operating points, followed by tne weignt removal 
transient run for tine domain validation of th2 linear 


heave-pitch model. 


It should be посед that further study of the 6 DOF 
model is needed to evaluate the Zesults of the two 
changes and to доге closely establish the arbitrary 
initial location of the center of pressure. The large 
amount of CPU and turn-around program times required to 
run tae nonlinear model precluded including that study 
in this work. Also, for tae nonlinear model, 
X-direction thrust was held constant during the run with 
the result that speed was 30.38 knots and still 
increasing at the end of the five second run. This 


effect will be further liscussed in the next section. 


C. COMPARISON OF RESULTS 


The graphical output for the five second cuns are 
shown іп figures 5 through 8 tor the linear nodel and 
figures 9 through 12 foc the nonlinear model. The fixed 
integration step size was 0.005 sec for the linear model 


and 0.001 sec for the nonlinear model. 


(see fig. 5 and 9| 
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In general both models exhibited a rapid initial 
reduction in pressure followed Бу а smooth transient 
back towards the final operating pressure. For the 
linear model the minimum occured at 100 milli-seconds, 
and 110 milli-seconds for the 5 DOF model, with the 
respective pressures being 21.38 and 21.30. Tae final 
values were 24.31 and 24.08 for the linear and nonlinear 


model respectively. 


AS сап be seen from the graphical data excellent 


agreement was acheived between the two models. 
2. Draft Transient 


(see fig. 6 and 10) 


Both models exhibit a smooth exponential type 
transient towards the final, lower steady-state draft. 
The linear model transient settles towards a lower value 
(6.04 inches as compared to 6.44) which is in part 
attributed то а difference in modeling of the seal 
forces. The nonlinear ı9del uses a complex conbination 
of hydrodynamic and aydrostatic lift where the linear 
model useS a Simple pressure differential effect. Also, 
Since the nonlinear model is settling to a significantly 
lower pitch angle (sse para. 4 below), the planing force 


effect is reduced in that model. 


ration 





= 
=a 


(see fig. 7 and 11) 


Both models show a rapid positive transient in 


acceleration that closely follows the negative transient 
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in plenum pressure. Both models reached the initial 
peak positive C.G. acceleration at approximately 20 
milli-seconds after the plenum pressure completed the 
negative transient. The linear model appears to have 
greater coupling between pitch and heave looking at the 
С.С. acceleration, but this is not actually supportable 
since the linear model has greater pitch amplitude 
excursions about the exporential pitch transient to the 
lower steady state value (see fig. 8 and 12.) ТП 
general the two models show excellent agreement іп 


transient and damped C.3. acceleration behavior. 


4. Pitch Transient 


(see fig. 8 and 12) 


It is in this area that the two models are in 
greatest disagreement. Both exhibit a general damped 


sinusoidal behavior about a center exponential transient 


to a lower steady state value. In both 23525 the 
initial excursion is downward.’ The natural pitch 
fequency of linear model is approximately 4.3 radians 


per second, whereas the nonlinear model exhibits a 
natural pitch frequency of 6.98 radians рес second. 
There is felt to be one major factor contributing to 
this Large discrepancy, which is tha inclusion of the 
planing force in the nonlinear model. Data presented by 
Thaler and Gerba 114] showed a natural pitch frequency of 
5.8 radians per second and much higher oscillations 
about the exponential transiant. It is felt that the 
Simplifiad planing forces of the linear 219421 аге 
included in the much more complex hydrodynamic equations 
of the nonlinear model and are therefore redundant in 
that model. Due to tine constraints, that avenue of 


investigation is left for future studies. The linear 
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model is tending towards a steady state value on the 
order or 0.45 degrees whereas the nonlinear model is 
tending towards approximately 0.30 degrees. This large 
discrepancy is partially due to the redundant planing 
force in the nonlinear model and the fact that speed is 
also increasing in the nonlinear model which tends to 
accentuate the planing force and force the craft to 
pitch down. 

In the next chapter a signal flow graph will be 
developed and Mason's Sain Rule applied to explain the 
linear model characteristic cesponse in a more 


analytical fashion. 
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А. INTRODUCTION 


The State Variables defined in Chapt IT are used to 
develop a signal flow graph of the craft heave-pitch 
dynamics. Mason's Gain rule is applied to the signal 
flow graph to derive analytical S-domain polynomials 
expressing the transfer functions for pitch moment into 
pitch angle, heave acceleration into draft and the two 
cross gains. Using typical gains calculated in the tine 
domain program developed in Chapt III the roots of these 


polynomials are located and it is shown that pitch 


dynamics may be very closely approximat2i by a 
second-order systen. Heave dynamics are also 
simplified. 


Using a simplified apdoroach to the frequency response 
of the craft to non calm-water conditions, the basics 
are developed for later sea-state model work. A dynamic 
system modeling program available at tae Naval 
Postgraduate School is used to generate Bode plots for 
the craft's frequency response in Heave and Pitch under 


the simplifying assumptions. 


В. SIGNAL FLOW GRAPH ANALYSIS 


Fig ШЕ was generated using the sensitivity 
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coefficients defined in Chapter III. From that graph the 
following feed-forward and feedback loops and paths were 
identified: 

1. Feed Forward, U(THETA) into THETA. 


P1 - (1/IYY)*(1/S**2) 


2. Feed Forward, U(Z) into Z. 


Р2 


(1/MASS)*(1/S**2) 


3. Feed Forward, U(T3ETA) into 2. 


P3 (1/IYY)*DZIH*(1/Sx**u) 


4. Feed Forward, U(Z) into THETA. 


рц (1/4ASS) *DTHZ*(1/S**4) 


5.Feedback Loops. 


11 = DTETH/S**2 

12 = DTHDTH/S 

L3 - (DZIH)*(DTHZ)/S**31 
L4 = (DZMB) *(DMBZ) /S **3 
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15 = DZZ/S**2 
L6 = DMBMB/S 
L7 = (DZTH) *(D4BZ) * (DTY4B) /S**5 


С. МАЅОМ* 5 GAIN RULE APPLIED 


Mason's gain rule for a closed multi-loop system is: 


Т = SUM(PK*DELTAK) /DELTA 


where, 

T = transmission gain from an input to an output 
node. 

DELTA = the graph determinant 

DELTAK = cofactor of the kth path 

Pk = gain of the direct path from input to output 


Using Mason's gain rale th? graph determinant (DELTA) 


рә tOund to be of the form: 


D (s) = (85*S**5 + Вужс̧жжц  B3J*s**3 c B2*s**2 + В1ж5 + 


ВО иж 


The forward direct path J (THETA) to THETA is: 


N1(S) = (1/LYY) * (1/3 **5) *(A3*S**3 + A2*S**2 + A1*S + AO) 


The forward path U(Z) to Z was found to be: 
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N2(s) = (1/MASS) *(1/S**5) *(A3*S**3 + A2*S**2 + A1*s + 


А0) 


The forward direct path U(THETA) to 2 15: 


N3(S) = (1/EYY) *DZTH*(S + AO) /S**4 


The forward path for draft coupling into 


to THETA is: 


NO = (1/MASS) *DTHZ* (S + AO) /S**4 


D. SOLUTION OF THE INDIVIDUAL POLYNOMIALS 


The gains used below were taken from a 


domain run as developed in chapter IV. 


DZZ = -2178.0 
DZTH = -121.3 
DZMB = -5499.0 
DTHZ - 28.78 
ОТНТН = -22.31 
DTHMB = 72.1 
БИН = -28.79 


DMBMB = -73.11 


1. Solution of th 


іш 


scaph Determinant 


pitch  (U(Z) 


typical tine 


Inserting the numerical values for the typical 


sensitivity coefficients listed above and using only the 
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dominant terms where possible, the roots of th2 graph 


determinant were found to be: 


DELTA = (S**2 + 71.8*S + 2066)*(S % 0.5086) «(5%%2 % 
0.72%5 %22.74) /5%ж5 


These factors represent a pair of complex poles 
located at -35.71+ 127.88 , а real pole located at 
-0.5085 , а pair of sonplex conjugat2 roots located at 
Е)- 350 % 2 4.755 , and à fifth order zero located at tne 


origin. 


25 оле о or As, 


№1 (5) = (1/1ҮҮ) * (5 + 0.427)« (5##2 + 72.68#5 + 2146) /5%%5 


This polynomial factors into a real zer» located 
at -0.427 , a pair of coaplex conjugate zeros located at 
510.38 + 3 28.74 , and a fifta order pole locatei at the 


origin. 


ВСИ сооп Of N2 (5) 


N2(S) = (1/MASS) *(S**2 * 0.7721*S * 22.35) *(S * 
72.38) /S**5 


Which factors into a pair of complex conjugate 
zeros located at -0.3856 + j 4.712 , a real zero at 


-72.38 , and a fifth order pole located at the origin. 
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N3(S) » (1/IYY) *DZTH*(S + 73.11) /5**% 


Which represents a real zero at -73.11, and a 


fourth order pole at the origin. 


5.  Soiution of N5 (S) 


— 


N4 (S) = (1/MASS) *DTHZ*(S + 73.11) /S**y 


Which has the same bole and zero locations as 
N3 (S). 


EAS EMPLIFICALION OF TRANSFER FUNCITIONS 


T1(S) = THETA(S) /U (THETA,S) 


Therefore from Mason's Gain Rule: 


T1(S) N1 (S) /DELTA 


T1(S) (iy са пиеш ее + О. 7 
(5 + 35.1 + ј 27.84) «(5 + 0.508) * (5 + 0.36 + ј 4.755) 


73 





To simplify let the complex conjugate zeros at 
-36.34 + j28.7 cancel the poles at -35.9 + 327.84 , and 
the real zero at -0.427 cancell the pole at -0.508. 


The result is the simplified seconá order systen: 
T1(S) = (1/TYY)/(5**2 ж 0.72*s + 22.74) 


Inspection of the time domain pitch response with 
respect to natural frequeny and damping confirms this 


basic underdamped second order response. 


Да тосу 


12 (5) 2 (5) /U (2,5) 


72 (8) (1/МА55) * (5 + 0.3856 + j4.712)*(S + 72.38) / 
(S**2 %35.1 * j 27.84) *(s * 0.5085) *((S +0.35 + 
34.755) 


Canceling the pole at -0.36 + j 4.755 wich the 
zero at -0.3856 + 34.712 results in the simplified 


response: 


T2(5) » (1/MASS)*(S * 72.33)/(S * 0.508) *(S**2 + 71.8%5 
* 2066) 


Further simplification can be made if (5 * 
ПІ е + 0.508) 25 3pprpoxsmated by 72.33/(S + 0.508) 
(which is valid for longer time period considerations.) 


Аз а result Т2 (5) becomes: 
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T2(S) 7 (72.33/MASS)/(8*0.508) *(S**2 « 71.8=5 + 2066) 


sn ceh раз poles located at -0,508, апі -35.9% 127,83 
This is in quite close agreement with the results 
obtained by Thaler and 3erba 14 foc the sama airflow 
condition (Ка + 1.0) for their heave-only linearized 
model. Their results were a real pole at -0.33 and a 


complex pair at -35.1+ j34.4 . This is of the form: 
T(S) = K/(S + 0.38) *(S**2 + 70.2%5 + 2415) 


ІН из (и 


T3(S) 7» Z(S) /U(THETÀ) 


Here it is noted that thera are no pole zero pairs 
that conveniently cancel so that tne transfer function 


becomes: 


T3(S) » (1/IYY)*DZTH*S*(S + 72.33) / 
(S + 35.1 + 5 27.84) *(s+ 0.5085) *(S + 0.36 + 34.775) 


Making the simplification used in ches previous 


section, the result iss 


Tomo = га ломитих) "ратни 5 /(5>7:2 в 71.985 + 2055) * 
(5 + 0.508) *(S**2 + 0.72*5 + 22.74) 


Which can be further Simplified by lookiag at the 
long term response vice initial few milli-seconds to 
O7 (S + О7508)) = 2 sọ tnat: 
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T3(S) -» (142/IYY)*DZTH/(s**2 * 71.8*s + 2056) *(5%ж2 + 
0.72*s * 22.74) 


à. TS) 


Since the transfer fundtions are the same except 
for the gains the solutions are the same except for the 


gains: 


T4 (S) = (1/MASS) *DTHZ * S*(S + 72.33) /DELTA 


where the same simplifying assumptions may be applied. 


Е. SIMPLIFIED BODE PLOT ANALYSIS 


1. Simplifying Assumption 


(see fig. 12) 


The assumption is made that the craft linearized 
frequency response transfer functions can be separated 


Ameo сео transfer functions. 


The first block E ен) expresses the 
disturbing moments ani accelerations generated by the 
crart interacting with the waves. This function was 
developed by Thaler and Gerba M4] for the heave-only 
model. In that Simpler model the function G1 was a 
complex function of encounter frequency, angle of 
incidence, craft length compared to encounter wavelength 
etc.. This portion of the overall frequency response of 


this model is left for future work due to time 
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eonstraints, 


The second block 32(We) represents the heave-pitch 
dynamics of the model as developed in chapters II, III, 
and IV, which respond to the disturbing moments and 
accelerations from block 61. 


The Transfer functions T1 and T2  (unsinplifizd) 
were Simulated using the IBM-360 simulation language DSL 
with S = jW with the following results. 


2. Bode Plots for Pitch Dynamics 


(see figs. 15 and 316) 


Figure 15 represents the normalized (numerator 
gain constants set to 1) gain in dB's versus Omega. 
Note that the abscissa scale is Ouega in powers of 10 
since this is a convenient means to achieve a 
logarithmic scale In the X-Y plotter output. As can be 
Seen the response strongly supports the simplification 
to second-order system. These plots represent the full 
transfer function Т1 (3м) vice the simplified systen. 
Figure 16 further supports the simplification as the 
phase plot starts at zero degrees and shows only a very 
slight positive increase prior to the rapid change 
toward -180 characteristic of an underdamped systen. 
Also note that above the damped natural frequency the 


magnitude plot falls at -12 dB per octave. 


о осе РП оез foc Heave Dynamics 


(see figs. 17 and 18) 
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Figure 17 represents the normalized na3nitude 
versus Omega for the ћгауе dynamics. The system is 
characteristic of an overdamped third order system with 
small gain increases at the craft's natural pitch 
frequency, which represents the pitch dynamics feeding 
back into heave dynamics. However the gain magnitude 
curve has already fallen 20 dB below it's lower 
frequency response gain before the pitch coupling begins 
to manifest itself. As can be seen in fig. 18 the phase 
plot shows the same third огаэг pole response with 
Slight second order zero influence at the natural pitch 


frequency. 
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А5 shown in Chapters IV andi V the simplified 
Pizch-Heave model developed herein is a very reasonable 
sinulation of the XR-3 test crait for the given 
operating conditions. Heave, “Pressure, and СО. 
Acceleration were in very close agreement. Pitch on the 
other hand showed some difference in natural frequency 
and damping factor. It is felt tnat this difference is 
primarily attributable to the inclusion of a redundant 
planing force calculation en the 6 DOF - model. 
Historically the results of the 6 DOF computations 14 
for pitch response have been very Similar to tha results 
for the linear model, in both natural frequency and 
damping characteristics. This can be attributed to the 
fact that the additional planing force did not 2xist in 


the original 6 DOF program. 


Since one of the major concerns with the 6 DOF model 
15 the excessive computational times involved in 
sea-state simulations the conclusion may be drawn that 
eventual substitution of the sidewall and seal force 
linear aodel eguation in the 5 DOF progran will 
substantially reduce 320 time, with some resultant loss 
іл сСоприсасізпаі accuracy. Ву way of comparison, 
disregarding compiling time and output time, the linear 
model reguired approximately 9 seconds of CPU time for 
the 5 second simulation while the 6 DOF model required 9 
Minutes of СРЈ for the 5 second run. This 15 а 


computation time improvenent on tne order of 60 times. 


From the time-domain plenum gauge pressure data it can 
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be seen that the dynamic variation in effective plenun 
roof area has a noticeable effect on the plenum pressure 
transient behavior as conpared to previous results 14. 
The initial pressure down transient was faster, dropped 
to a lower value , and began it's recovery more quickly. 
Therefore the plenum roof dynamic variations strongly 
affects the C.G. acceleration characteristics of the 
craft and therefore it's transient behavior. 

The dynamic center of pressure variation developed 
herein proved to be a very draft sensitive ph2nomenon. 
Although the length of the other lever arms as nodeled 
herein change by a similar amount, percentage of 
variation is radically different. For the seals the 
variation is on the order of inches in 10 fest, where 
the variation in lever arm for plenum pressur2 lift is 
on the order of inches in a single foot (or less.) This 
implies that dynamic C.P, variation 15 th? lominant 
ractor in steady state pitch angle versus steady state 


dratt. 


Planing forces as included in this model are now in a 
form where the basic Macros ас rir ase cratt 
construction characteristics are readily recognized. 
This force provides part of the answer to craft 
pitch-down characteristics as initially investijated by 
Reidel [S], with the remainder of the moment being 
developed by the С.Р. variation for pressure. The 
planing force provides a clue as to how, during 
construction, the craft's pitch amplitude zxzursion 
characteristics may be controlled. 


Several areas for future investigation have become 


apparent in this work. 


1. Verification of the effect of dynamic center of 
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pressure variation and planing forces as incorporated in 


the 6 DOF model over the entire speed operating range. 


2. Incorporate the simple modeling for sidewalls and 
seals into the 6 DOF model and investigate the trade-off 
between CPU time and conputational accuracy. 


3. A closer look at the seal forces used herein. Xi 
is possible that a simple gain factor needs to be 
included in these equations in order to optimize the 
Steady-state seal forces with respect to actual load 
data as measured during experimental test runs of the 
craft 10. 


4. Incorporate sea state forces and moments into the 
Simplified model by simulation of the forces and monents 
that are generated by craft interaction with the sea. 
Evaluate the accuracy of the method used by comparison 
With data generated by the simple model and th2 6 DOF 


model in sea-state conditions. 


5. Analytical investigation of the craft physical 
characteristics that have predominant influence on the 
natural pitch frequency and damping by using the gain 
terms from the results of the signal flow graph analysis 


done in Chapter V. 


о. Use the simpiifiel model to develop a ride-conmfort 
control system. Test the resulting design on the 6 DOF 
model, and compare the results obtained with those 
suggested by RIAR [7]. It has been known ampirically 
for some time that the craft exhinited simple second 
order behavior in the pitch mode. Now that this fact 
has been analytically verified the implications for an 
automatic control system for pitch response аге 


powerful. A simple control system can be visualized 
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that does not involve zither complex optimal control 
system design work or hardware. A simple acceleration 
and rate feedback control system is implied, with the 
possibility of adding a speed input to the system to 
compensate for the localized linearity of the craft 
characteristics. 
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АРРЕНОТХ А 


LINEAR SES MODEL 958815 GUIDE 


AS General Comments 


The program as listed following Appendix B is intended 


for useage via the hot card reader at NPS. 


It is intended to run in Fortran С with Plotting 
option (PROC = FORTCLGP.) 


The program 15 structured on the assurption that the 
initial and final state variables аге іп equilibriun 
prior to data entry, but will run regardless. Га force 
and moment residual listings will indicate the value of 


imsalance, if any. 


The craft physical dimensions are alterable by 
changing the cards in the craft dimension section of the 
Main program, which are keyed to the craft geometry 
herein. 

В. Data Deck Setup. 


The /* and //GO.SYSIN DD * cards are inserted 
following the end card in Subroutine VERSAP, followed by 
the data deck as described 


Peer Opteron сагп (ТӘрОТ), тогпас ІТ. ТРГОГ = 1 = 


ESSSATECCSX-Y plotter output. 
IPLOT = 0 = PRT PLOT output 
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Ж Ілес option carda (12157), Fornat I1 

ILIST = 1 = Forces and moments,with residuals and all 
sensitivity coefficients listed. 

ILIST = 0 = Input craft operating conditions only 
listed. 


3. Time and Integration variables. 

ТГ, ТР, NPIS, NSTEP Pormat (2F10.4,2110) 

Ti = problem start time (normally 0.0) 

TF = Problem stop time (seconds) 

NPTS = number of tabular output points desired. 

NSTEP = number of integration steps. 

note: it is best to make NSTEP an integer multiple 
of NPTS, and never less than NPTS. 

4.Planing coefficient and arbitrary center of pressure 
input, i.e. READ( ) AKP, PLCOEFF 

Format = 2F10.4 


5% Initial Conditions card. 

The initial craft conditions are read in in the 
following order: Weight (Lbs.), Draft (in), Pitch angle 
(Deg.)), EBBAR (ps£f), and speed (ips.) Format is 5F10.4 

6. Final Conditions. 

The final conditions are read in the following order 
and format: waight (LDS); Drift (їй), Pitch angle 
(Deg), and PBBAR(psf), Гогпае 4F10.4. Since this is a 
constant speed model, initial and final speed are 
assumed 2qual. 


Је Label Cards. 


Label cards for VERSATEC output are the Last 10 cards. 
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There are five graphs in this mode. . Draft(in), C.G. 
Acceleration (G's), Pitch Angle (Deg.), Pitch 
Acceleration (Deg/SEI**2), and PBBAR (Each graph 
requires two label cards, Format (6A8), which print on 


two lines below the graph. The graphs are output in the 
order listed above. 
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APPENDIX B 


DRAWINGS AND FIGURES 
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